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Abstract 

We report a precise measurement of the lifetime. The data were taken by the 
SELEX (E781) experiment using 600 GeV/c £~, ir~ andp beams. The measurement 
has been made using 83 ± 19 reconstructed Q,® in the Q~ tt~ tt + tt + and f2~7r + decay 
modes. The lifetime of the is measured to be 65 ± 13(stat) ± 9(sys) fs. 
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1 Introduction 



Several experiments [1-6] in the last years have detected the ground state 
as well as an excited state. Recently at Fermilab the photoproduction experi- 
ment, FOCUS, reported an observation of a sample of 64 Q° events and they 
measured its lifetime as 72 ± 13 ± 11 fs [7]. The experiment WA89 published 
an lifetime measurement of 55^ \l j^i & from a sample of 86 events [8] . The 
present world average is 69 ± 12 fs as reported in Ref. [9]. Clearly additional 
measurements of lifetime as well as the branching ratios with more statistical 
accuracy are needed to test theoretical models [10,11]. 

In this letter we report the results of a new measurement of the lifetime based 
on data from the hadropro duct ion experiment SELEX (E781) at Fermilab. 
The measurement is based on a sample of 83 ± 19 fully reconstructed from 
15.3 x 10 9 hadronic interactions. 



The SELEX detector at Fermilab is a 3-stage magnetic spectrometer. The 
negatively charged 600GeV/c beam contains nearly equal fractions of E and 
7i. The positive beam contains 92 % protons. Beam particles are identified 
as a baryon or a pion by a transition radiation detector. The spectrome- 
ter was designed to study charm production in the forward hemisphere with 
good mass and decay vertex resolution for charm momenta in the range 
100 — 500 GeV/c. Five interaction targets (2 Cu and 3 C) have a total target 
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thickness of 4.2%\i nt for protons. The targets are spaced by 1.5cm. Down- 
stream of the targets there are 20 Silicon Strip Detectors (SSD) with a strip 
pitch of 20 — 25 /im oriented in X, Y, U and V views. The first spectrom- 
eter level has three MultiWire Proportional Chambers (MWPC) with 3 mm 
wire spacing and 2 x 2 m 2 area downstream of bending magnet Ml. The sec- 
ond spectrometer level has 7 MWPCs with 2 mm wire spacing downstream 
of the second bending magnet M2. Each chamber has two sensitive planes 
in two orthogonal projections. The scattered-particle spectrometers have mo- 
mentum cutoffs of 2.5GeV/c and 15GeV/c respectively. Typical momentum 
resolution for a 100 GeV/c track is 0.5 %. A Ring-Imaging Cherenkov detector 
(RICH) [12], filled with Neon at room temperature and pressure, provides sin- 
gle track ring radius resolution of 1.4% and 2 a K/n separation up to about 
165GeV/c. A layout of the spectrometer can be found elsewhere [13]. 



2 Reconstruction of hyperons f2, S and £ 

The decays studied here have a hyperon in the final state. Hyperons that 
decay upstream of or within the Ml magnet are called Kink tracks. They are 
characterized by one charged track that decays to another charged particle 
and a neutral particle undetected by the spectrometer. Such Kink tracks dif- 
fer from the majority of spectrometer tracks in that the vertex silicon track 
segment does not link to straight line tracks segments measured in the spec- 
trometers after Ml and/or M2. The Kink reconstruction algorithm examines 
all unlinked Vertex SSD track segments that point to the Ml magnet aperture 
and tries to match each unlinked segment with an unlinked downstream track 
measured in the M1/M2 spectrometer, using momentum-energy conservation 
with the hypothesis of a specific hyperon decay. The momentum of the parent 
hyperon (f^, S 1 * 1 or S ± ) is calculated using the assumed decay hypothesis. 
The daughter K~ in Q~ decays must be RICH identified with the likelihood 
to be a K exceeding that of its being a it. 



3 Data set and charm selection 

The charm trigger is very loose. It requires a valid beam track, at least 

4 charged secondaries in the forward 150 mrad cone, and two hodoscope hits 
after the second bending magnet from tracks of charge opposite to that of the 
beam. We triggered on about 1/3 of all inelastic interactions. A computational 
filter linked MWPC tracks having momenta > 15 GeV/c to hits in the vertex 
silicon and made a full reconstruction of these tracks together with a beam 
track to form primary and secondary vertices in the event. Events consistent 
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with only a primary vertex are not saved. About 1/8 of all triggers are written 
to tape, for a final sample of about 10 9 events. 

In the full analysis the vertex reconstruction was repeated with tracks of all 
momenta. The RICH detector identified charged tracks above 25GeV/c. Re- 
sults reported here come from a second pass reconstruction through the data, 
using a production code optimized for hyperon reconstruction. 

To separate the signal from the non-charm background we require that: (i) the 
spatial separation L between the reconstructed production and decay vertices 
exceeds 6 times the combined error a^, (ii) each decay track, extrapolated 
to the primary vertex z position, must miss by a transverse distance s > 2.5 
times its error a s , (hi) each candidate hyperon track, extrapolated to the Kink 
vertex z position, must have a good vertex quality (x 2 /NDOF < 5), (iv) the 
secondary vertex must lie outside any target material by at least 0.05 cm, and 
(v) decays must occur within a fiducial region. 

The total transverse momentum of pions from the fi~7r + 7r + 7r~ decay mode 
must be greater than 0.35GeV/c with respect to the direction. This cut 
optimizes the signal to background ratio. We require a minimum n momen- 
tum of 8GeV/c to reduce the number of fake invariant mass combinations. 
There are 107 ±22 candidates in three decay channels: f2~7r + 7r + 7r~, Q~n + , 
and E~K~it + it + . Details of the Q° mass measurement will be reported else- 
where [14]. In the E~K~it + it + mode the signal is small and the signal to noise 
ratio is poor. We've chosen not to include it in the lifetime measurement. The 
invariant mass distributions with tighter cuts for lifetime measurement are 
shown in fig. 1 for the two decay modes used here. 



4 Lifetime evaluation using a maximum likelihood fit 

The average combined error gl on the primary and secondary vertices and the 
average momentum give a proper time resolution of 16 fs, about 23 % of 
tqo. We used a binned maximum likelihood fitting technique with 5fs width 
bins to determine the lifetime. The fit was applied to a reduced proper 
time distribution, t* = M(L — L min )/pc where M is the reconstructed charm 
mass, p the reconstructed momentum, L the measured vertex separation and 
L m i n the minimum L for each event to pass all the imposed selection cuts. 
L min is determined event-by-event, along with the acceptance, by the proce- 
dure described below. We fitted all events with t* < 600 fs in the mass range 
2.685 GeV/c 2 < M(Q° C ) < 2.725 GeV/c 2 , ± 2.5 a from the tt° c central mass 
value. 

The probability density function is: 
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Fig. 1. Invariant mass distribution for (a) 0^? — ^ 7r + 7r 7r + , (b) — > f2 tt + . The 
shaded regions show the signal and sideband regions. 

f(r n o,T B l,T B 2,fB,fc,t*) = 



(1 - f B )N s 



e{t*)r n o 



+ f B N s B(t*) 



where 



P -t*/r B1 p-t*/r B2 
B(f) = fc + (l-/c)- 



TB1 



(2) 



The five parameters are: tqo (Q° lifetime), r B1 , r B2 (background lifetimes), 
/b (background fraction in the signal region) and f c (background splitting 
function). Ns is the total number of events in the signal region. 

The function is the sum of a term for the exponential decay corrected 
by the acceptance function e(t*) plus a background function B(t*) consist- 
ing of two exponentials to describe the strong decays and charm decays re- 
spectively. Its parameters were determined from the t* distribution from the 
sidebands. Together the mass widths of the sideband background win- 
dows, 2.610 GeV/c 2 < M(Q° C ) < 2.625 GeV/c 2 and 2.760 GeV/c 2 < M(Q° C ) < 
2.785 GeV/c 2 was equal to the signal mass window. 



The proper-time-dependent acceptance e(t*) is independent of spectrometer 
features after the first magnet, e.g., RICH efficiency and tracking efficiency 
These efficiencies affect only the overall number of events detected. The proper 
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Table 1 

Lifetime fit results for the two 0° decay modes analyzed. The errors are only sta- 
tistical. The signal yields from the fits to the mass plots in Fig. 2 are also shown. 



Fit results 


7T + 7T + 7r 


Q-TT+ 


T n» [fs] 


62.6 ± 22.0 


65.8 ± 16.0 


TBI N 


15.6 ±6.2 


10.1 ±3.3 


TB2 M 


388.2 ± 27.0 


281.4 ± 22.5 


Signal 


34 ± 12 


23 ±9 


background 


84 ± 13 


81 ±6 


Signal yield 


47 ± 16 


36 ±11 



time distribution of these events depends crucially on vertex reconstruction. 
To evaluate e(t*) we generated events with a (1 — xf) 3 distribution and 
decayed them using the QQ package [15]. We embedded these generated decays 
into real data events and reconstructed the embedded decays with the offline 
package including multiple Coulomb scattering in the spectrometer and the 
measured detector performance. The correction function was evaluated as the 
fraction of the embedded events passing the selection cuts. 

Figure 2 shows the overall fits to the data distributions as a function of re- 
duced proper time for f2~7r + 7r~7r + and decay modes. It also shows the 
acceptance function e(t*), which does not differ significantly from unity and is 
constant. This is due to the fact L min is chosen 6 times the combined error o^. 
With SELEX's very high momentum and excellent resolution this cut removes 
only the first ~ 1.5 lifetimes from the sample. 

Table 1 summarizes the lifetime fit results and the signal yields. We measure 
an average lifetime 65 ± 13 fs. The uncertainties are statistical only, evaluated 
where — ln£ increases by 0.5. 



5 Systematic errors 

The systematic uncertainties for the lifetime analysis are listed in Table 2 
and described below. Lifetime shifts due to reconstruction errors have been 
well studied in our D° and A c work, with an order of magnitude higher statis- 
tics [16,17]. Because of the high redundancy and good precision of the sili- 
con vertex detector, vertex mismeasurement effects are small at all momenta. 
Proper time assignment depends on correct momentum determination. The 
SELEX momentum error is less than 0.5 % in all cases. We assign a maximum 
systematic error from proper time mismeasurement of 1 fs. The acceptance 
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Fig. 2. Corrected reduced proper time distribution for events in the window 
2706 ±20 MeV/c 2 (full boxes) and the results from the maximum likelihood fit (solid 
curve) for a) f2 - 7r + 7r + 7r~; b) 7r + . The dashed curve shows the 0° proper lifetime, 
the dashed-dot curve the fitted background and the dotted curve the acceptance. 

Table 2 

Systematic error contributions in fs for the two decay modes analyzed. 



Source of uncertainty 


Q 7r + 7T 7T + , Q 7T + 


Vertex reconstruction 


< lfs 


Acceptance function 


1.75 fs 


Fit procedure 


8.5 fs 


Total systematic error 


8.6 fs 



function used in the fit was parameterized with a 1st and 2nd order polyno- 
mial. The difference in lifetime result is 1.75 fs. No significant difference in the 
lifetime correction function was found when we changed the n value of the xp 
distribution from 3 to 1. We varied the width of the sidebands and the bin 
size independently. The systematic error due to the fit procedure is 8.5 fs. 



6 Conclusions 



We have made a new measurement of the VP C lifetime in two independent 
decay channels, f2~7r + 7r + 7r~, fi~7r + , using a maximum likelihood fit. SELEX 
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measures the f2° lifetime to be tqo = 65 ± 13(stat) ± 9(sys) fs. Our result is 
in excellent agreement with the world average [9] and with the recent results 
published by the FOCUS collaboration [7]. 

Acknowledgment 

The authors are indebted to the staff of Fermi National Accelerator Lab- 
oratory and for invaluable technical support from the staffs of collaborating 
institutions. This project was supported in part by Bundesministerium fur Bil- 
dung, Wissenschaft, Forschung und Technologie, Consejo Nacional de Ciencia 
y Tecnologia (CONACyT), Conselho Nacional de Desenvolvimento Cientifico 
e Tecnologico, Fondo de Apoyo a la Investigation (UASLP), Fundacao de 
Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), the Secretaria de 
Education Piiblica (Mexico) (grant number 2003-24-001-026), the Israel Sci- 
ence Foundation founded by the Israel Academy of Sciences and Humani- 
ties, Istituto Nazionale di Fisica Nucleare (INFN), the International Science 
Foundation (ISF), the National Science Foundation (Phy #9602178), NATO 
(grant CR6. 941058-1360/94), the Russian Academy of Science, the Russian 
Ministry of Science and Technology, the Russian Foundation for Basic Re- 
search: RFBR grant 05-02-17869, the Turkish Scientific and Technological 
Research Board (TUBITAK), the U.S. Department of Energy (DOE grant 
DE-FG02-91ER40664 and DOE contract number DE-AC02-76CHO3000), and 
the U.S.-Israel Binational Science Foundation (BSF). 

References 

[1] J. Stiewe, AIP Conf. Proc. 272 (1993) 1076. 

[2] P. L. Frabetti et al. [E687 Collaboration], Phys. Lett. B 338 (1994) 106. 

[3] P. L. Frabetti et al. [E687 Collaboration], Phys. Lett. B 300 (1993) 190. 

[4] D. Cronin-Hennessy et al. [CLEO Collaboration], Phys. Rev. Lett. 86 (2001) 
3730 [arXiv:hep-ex/0010035]. 

[5] P. L. Frabetti et al. [E687 Collaboration], Phys. Lett. B 357 (1995) 678. 

[6] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 97 (2006) 232001 
[arXiv:hep-ex/0608055]. 

[7] J. M. Link et al. [FOCUS Collaboration], Phys. Lett. B 561 (2003) 41 
[arXiv:hep-ex/0302033]. 

[8] M. I. Adamovich et al. [WA89 Collaboration], Phys. Lett. B 358 (1995) 151 
[arXiv:hep-ex/9507004]. 



9 



[9] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[10] S. Bianco, F. L. Fabbri, D. Benson and I. Bigi, Riv. Nuovo Cim. 26N7 (2003) 
1 [arXiv:hep-ex/0309021]. 

[11] H. Y. Cheng, Phys. Lett. B 289 (1992) 455. 

[12] J. Engelfried et al. [SELEX Collaboration], Nucl. Instrum. Meth. A 431 (1999) 
53 [arXiv:hep-ex/9811001]. 

[13] J.S. Russ et al. [SELEX Collaboration], in Proceedings of the 29th International 
Conference on High Energy Physics, 1998, edited by A. Astbury et al. (Word 
Scientific, Singapore, 1998) Vol. II, p. 1259 [arXiv:hep-ex/9812031]. 

[14] SELEX Collaboration, to be published. 

[15] http : //www. Ins . Cornell . edu/public/CLEO/ sof t/QQ/ 

[16] AY. Kushnirenko, Ph.D. Thesis, Carnegie Mellon University, 2000 
(unpublished). FERMILAB-THESIS-2000-09. 

[17] A. Kushnirenko et al. [SELEX Collaboration], Phys. Rev. Lett. 86 (2001) 5243 
[arXiv:hep-ex/0010014] . 



10 



